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1 INTRODUCTION 



ABSTRACT 

We study a composition of normal and exotic matter which is required for a flat 
Emergent Universe scenario permitted by the equation of state (EOS)(p = Ap — Bp^) 
and predict the range of the permissible values for the parameters A and B to explore 
a physically viable cosmological model. The permitted values of the parameters are 
determined taking into account the H{z) — z data obtained from observations, a model 
independent BAO peak parameter and CMB shift parameter (WMAP7 data). It is 
found that although A can be very close to zero, most of the observations favours a 
small and negative A. As a consequence, the effective Equation of State parameter for 
this class of Emergent Universe solutions remains negative always. We also compared 
the magnitude {p{z)) vrs. redshift(z) curve obtained in the model with that obtained 
from the union compilation data. According to our analysis the class of Emergent 
Universe solutions considered here is not ruled out by the observations. 
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It is well known from the recent observations that the 
standard Big Bang model of cosmology fails to de- 
scribe the present accelerating phase of the universe. The 
model is also pleagued by a time like singularity in 
the past. Accelerating phase of the universe can, how- 
ever, be incorporated in a number of ways. A number 
of m odels e.g., models with modified theory of gravity 
i riou 20071), models with u nusual matt ers like Chaplygin 
gas I Bento. Bertolami fc Senll2002l : iBilic. Tupper fc Viollieii 
I2OO1I ). scalar and tachvon fields l|Lvthl 1200 J T are taken 
into account to accommodate present phase of acceleration. 
There are other models based on mostly non-equilibrium 
thermodynamics and Boltzman n formulatio n which do not 
require dark energy ( Zimdahl et al.l 1200 ll : iBalakin et al.l 
I2OO3I : Ibima. Silva fc Santosi. 20081 ). Only very recently some 
other models appeared in the literature which discusses 
cold dark matter (CDM) and CDM interactions as alter- 
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native to the ACDM mo del l|Lima. Jesus fc OliveirallioogI : 
iBasilakos fc PlionisI |2009| ). However, exploring singularity 
free cosmological models is an interesting area in cosmol- 
ogy and Emergent Universe scenario (EU) is one of the 
well known choices. A number of literature appeared which 
discussed EU model as it was free from initial singular- 
ity and the size of the universe was large enough so that 
quantum gravity effect s were not important (|Harriso 3 ll967l : 
lEUis fc Maartenslliool ). These models evolve from a static 
phase in the infinite past into an inflationary phase. The idea 
is in conformity with Lemaitre-Eddington concepts from 
early days of modern cosmology. If developed in a consis- 
tent manner an emergent universe model is capable of solv- 
ing some of the well known conceptual problems not under- 
stood in the Big-Bang model. A model of an ever-existing 
universe, which eventually enters into the standard Big Bang 
epoch at some stage and consistent with features known 
to us today is worth considering. Recently an inte r estin; 
EU model has been proposed by iMukheriee et al] (|200' 
which requires some exotic matter in addition to normal 
matter as cosmic fluid. The model has been explored in a 
flat universe as such universe is supported by recent obser- 
vations. Subsequently the EU model was taken up to ex- 
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amine the suita bility of implementing it in the context of 
various theories llBaneriee. Bandvopadhvav fc Chakrabort 



Dcbnath"2008l: iBeesham. Chervon fc Maharaill 200' 
The EOS needed for the model pro- 
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posed by IMukheriee et al] (|2006l ') is given by 



i 



Table 1. H{z)vs.z data 



Ap- Bp^ , 



(1) 



where A and B are unknown parameters of the theory. We 
note that different values of A and B pick up different com- 
position of matters which may lead to an EU model. A 
similar EOS was considered in the literature as a double 
component dark energy model (|Fabris et al.l [2007) where 
the model parameters are constrained from Type la super- 
nova data. The EOS considered by them is basically a spe- 
cial form of a more general EOS, p = Ap — Bp"; which 
ma y be seen as a realization of Chaplygin gas (with a < 
0) llBento. Bertolami fc SenI |2002| : iBilic. Tupper fc VioUieil 
l200lf r It may be mentioned here that Chaplygin gas models 
were introduced in cosmology as an interpolation between a 
matter dominated era and a de Sitter phas e. Later a modi - 
fied model of Chaplygin gas was proposed (|Liu et a l. 2005') 
to describe cosmological evolution. For example models like 
Modified Chaplygin gas serves well as an interpolat i on be - 
tween radiative era and ACDM era. iFabris et al] (|2007l ) 
showed that even with a > such interp olation is permis- 
sible and an EOS like one considered by IMukheriee et al] 
(|2006l ) may be considered as a phenomenological realiza- 
tion of string specific c onfiguration. The model proposed by 
IMukheriee et alj (|2006l ) developed the EU scenario in a very 
general way. Once one considers an EOS given by eq. (1), a 
class of EU solutions is permitted for B > 0. The authors 
however showed that the above EU solution are not permit- 
ted with a minimally coupled scalar field. Also, it was found 
that the EU scenario automatically admitted a composition 
of three kinds of matter energy density in the universe all 
having their own way of evolution. This is certainly an inter- 
esting issue keeping in mind that the model has the provision 
for a large class of possible dark energy and dark matter can- 
didate. It is thus worth to investigate the viability of such 
an EU model with the recent observational data. Neverthe- 
less we intend to explore the allowed range of values of the 
parameter A, for B > for a viable cosmological scenario 
permitted by observations. To determine the range of val- 
ues for A and B for a viable cosmological model permitted 
by observations, we adopt here two independent techniques: 

i i) Applying m inimization technique on H{z) vs. z data 
Stern et ahibood ). Here we use 9 data points given in Table 
1, (ii) using joint analysis of H{z) vs. z data and a model 
independent BAO peak parameter and (iii) using joint anal- 
ysis of H{z) vs. z data, BAO peak parameter and CMB 
shift parameter. We explore here the suitability of the model 
with the help of supernovae data (union compilation data) 
also. The plan of the paper is as follows : In sec. 2 relevant 
field equations are obtained from Einstein field equation. In 
sec. 3 the constraints on model parameters are determined 
from H{z) vs. z data. Subsequently in sec. 4 and sec. 5 we 
obtain more stringent constraints on model parameters in 
accordence with the joint analysis with model independent 
BAO peak parameter and CMB shift parameter. In sec. 5, 
we draw the ^(z) — z curve for our model to compare with 
that drawn using union compilation data (|Kowalaski et al] 
I2OO8I ). Finally, in sec. 6 we summarize our results and briefiy 
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Figure 1. (Colour Online)Constraints from H{z) vs. z for K = 
0.0100. (Solid) 99% and 95% (Dashing) contours. The best fit 
point is shown (0.0122, -0.0823). 



discuss the results. 



2 FIELD EQUATIONS 

We consider Friedmann-Robertson-Walker(FRW) metric 
(c=l), given by : 



dr^ 



1 — fcr^ 



- r^{d9^ + sin? 6 1 



(2) 



where fc = 0, -|-1(— 1) is the curvature parameter in the spa- 
tial section representing flat, closed (open) universe and a{t) 
is the scale factor of the universe and r, 9, are the di- 
mensionless comoving co-ordinates. The energy conservation 
equation is given by 

J+3H(p-t-p) = 0, (3) 

where p, p and H axe respectively pressure, energy density, 
Hubble parameter. We express Hubble parameter {H) in 
terms of redshift parameter [z) which is given by 
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1 dz 



(4) 



1 + z dt 

Since components of matter (baryon) and Dark energy (ex- 
otic matter) are conserved separately, we use energy conser- 
vation equation together with EOS given by (1) to determine 
the expression for the energy density. Consequently eq. (3) 
yields: 



pemu — 



B 



1 + A 



+ 



K 



(5) 



where K is an integration constant which is required to be 
positive definite. From eq. (5) it is evident that the energy 
density is composed of three different terms, where a con- 
stant term {{j^)"^) may be identified with a cosmological 
constant and the other two terms are identified with two 
kinds of fluids determined by the parameters A. For sim- 
plicity, eq. (5) can be rewritten as. 



Pemu — PerUQ 



As + 



As 



,3{A-|-l)/2 



(6) 



where As 



l+A 



and 



K 
A+l 



' Pemo 



B 
A+l- 



Using the Friedmann equation we express H in terms of 
redshift parameter (2) for the model, which is given by 



Hiz) 



(l-n(>o)[ 



B + K ^ ^ ' ^ ' 

— f , where Q is composed of baryon 
and exotic fluids. Qbg represents baryon energy density and 
f2emo repre sents the exotic fluid densi ty. Here we consider 
Qba = 0.04 l|Dev. Alcainz fc Jain|[2007l ). 



with fl — Qtg + f2 



3 H(Z)-Z DATA AS A CONSTRAINING TOOL 

The EU considered here is implemented in a flat universe, 
consequently we consider baryonic matter and the exotic 
matter in a fiat Friedmann universe to the allowed range 
of the parameters. The Hubble parameter given by (7) is 
a function of a number of variables, consequently we can 
re-write eq. (7) as : 



H^{Ho,A,B,K,z) = HoE^{A,B,K,z), 
where 

E{A,B,K,z) = [ftt„{l + zf + 

.B + K{1^ 



3(A + 1) 



B + K 



(8) 



(9) 



For a given value of K, the best fit values for the unknown 
parameters of the model, namely A and B are determined 



by minimizing a Xh- 



X%-.{Ho,A,B,K,z) = J2 



function which is given below. 

{H{Ho,A,B,K,z)-Hots{z)f 



(10) 



where Hohs(z) is the observed Hubble parameter at redshift 
z and Oz is the error associated with that particular ob- 
servation. Since we are interested in determining the model 
parameters, Hq is not important for our analysis. So we 
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Figure 2. (Colour Online) Constraints from joint analysis with 
H{z)-z data BAG peak parameter for K = 0.0101. 99% (Solid) 
and 95% (Dashing) confidence level are shown in the figure along 
with the best fit value(0.0094,-0.1573) 



marginalize over Ho to get the probability distribution func- 
tion in terms oi A,B,K only, which is given by 

'xl-.{Ho,A,B,K,z 



L{A,B,K)^ dHoPiHo)exp 



where P{Ho) is the prior distribution function for the 
present Hubble constant. Here we consider Gaussian pri- 
ors, Ho = 72 ± 8 . One can minimize x^ by maximizing the 
function L{A, B , K). We fix K at the best fitted value and 
contours in A-B plane are drawn. However, fixing of K is al- 
lowed as we are interested to obtain range of A and B which 
is related to the EOS given by eq. {!)■ K enters in the theory 
as an integration constant which is always positive {K > 0). 
In fig. 1 we draw 99% and 95% contours on A-B plane. We 
see that within 99% confidence -0.5949 A f .663 and 
-0.0022 ^ B s; 0.3189. Of course theoretically we must 
have B ^ to obtain an EU solution for the model we are 
considering here. 



4 JOINT ANALYSIS OF H{Z)-Z WITH BAO 
PEAK PARAMETER 

In this section we use the technique adopted by 
lEiscnstei n et al.l (|2005l l to explore the parameter A (which is 
determined by A). The value of A is independent of the cos- 
mological models, which for a fiat universe can be expressed 



(11) 



A- 



Jo 



^1 dz 
E(z) 



2/3 



£(zi)l/3 



(12) 



where = + {1 - ^^(l " B/{K + B)f and 21 = 0.35 
and A = 0.469 ± 0.017. We define xIao = '"In °;t?/ , and 
for joint analysis we consider Xjoint = ^ 

The above joint analysis scheme with BAO sets new 
constraints on A and B, again which are upto 95% confi- 
dence level -0.3053 ^ 4 < -0.0306 and 0.0077 < B ^ 



(0.017)^ 

+ Xbao- 
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Figure 3. (Colour Online) Constraints from joint analysis of 
H{z)-z data, BAO peak parameter and CMB shift parameter for 
K = 0.0101. 99% (Solid) and 95% (Dashing) confidence level are 
shown in the figure along with the best fit value(0. 0102, -0.0176) 




Figure 4. (Colour Online) fi{z) vs. z curve comparism with su- 
pcrnovae data 



0.0116 and upto 99% confidence level -0.2757 ^ A 
-0.0500 and 0.0078 s: B ^ 0.0114. 



5 JOINT ANALYSIS WITH H{Z)-Z , BAO 
PEAK PARAMETER AND CMB SHIFT 
PARAMETER (7^) 

CMB shift parameter (TZ) is given by 



7^ = 




(13) 



z at the surface of last scattering. The 
gives us 7^ = 1.726 ± 0.018 at 2 = 



where zi 

WMAP 7 data 

1091.3 d Komatsu et~al]|2010l ). Thus we consider Xcmb 



(T;,-1.726) 



(0.018) 



^, with XTo 



xl- 



+ X%AO + Xcmb which im- 
pose additional constraints on the model parameters. The 
statistical analysis with XTot further tightens up the bounds 
on A and B. In fig. 3, 95% and 99% contours are plotted 
on A-B plane. We determine constraints from this analy- 
sis: within 95% confidence level -0.037 ^ A 0.008 and 
0.008 ^ B < 0.013. However, within 99% confidence limit 
we get -0.034 < A ^ -0.0114 and 0.012 sC B 0.009. The 
best fit value obtained here is given hy A = —0.0219 and 
B = 0.0103. Finally we draw a magnitude {n{z)) vs. red- 
shift {z) curve for our model with the best fit values of A, 
B and K and also show th e same curve drawn fr om union 
compilation data for SNela (|Kowalaski et al.]|2008l ) in fig. 4. 



6 DISCUSSIONS 

Theoretically the class of EU solutions considered here can 
be realized for a comp osition of difi'erent kinds of matter 
l|Mukheriee et al.l I2OO6I ) depending on the model parame- 
ters A and B. In this paper we determine allowed ranges 
for the model parameters (particularly those involved with 
the EOS i.e., A and B). Considering a EOS required for 
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Figure 5. (Colour Online) Plot of density parameter (O) for 
effective dark energy and effective matter content of Universe. 



a flat emergent universe we determine the constraints on 
the parameters using data available from cosmological ob- 
servations. We note from the analysis that the EOS that 
permits a class of EU solution, considered here, should con- 
tain exotic matter {A < Q, B > Q). This is certainly not 
ruled out by the theory itself. It may be pointed out here 
that H{z)-z data puts a bound on the model parameters 
which is further investigated in the light of the other obser- 
vational data such as BAO peak parameter value and CMB 
shift parameter. In the above we estimate model parame- 
ters {A, B and K). Most important point to be noted here 
that the later observations do not permit a positive value 
for the parameter A. Only small negative values seem to be 
allowed. Although positive A values are permitted when we 
consider H{z)-z data only but the best flt value is found 
to be negative. However the possibility that A ^ can not 
be entirely ruled out since our analysis permits values of A 
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Figure 6. (Colour Online) Effective EOS parameter for EU ( 
i^eff ) is plotted with redshift. Plot with the best fitted values 
of model parameters k = 0.0102, A = -0.0176 and B = 0.0103 
(Solid). Plot with values within 95% confidence (Dashing) and 
within 99% (Dotted). 



which are even very close to zero {A — 0) and the model 
may be realized in the presence of dust and dark energy. We 
also study the evolution of various cosmological parameters 
of the model. For example density parameter is on such im- 
portant parameter. We plot density parameter for effective 
dark energy and effective matter content of the universe with 
the redshift in fig. 5. We note that almost 80% of the present 
matter-energy content is effective dark energy and baryonic 
and nonbaryonic matter constitutes the remaining part. So, 
as far as present budget of Dark Energy and Dark matter is 
concerned, EU differes very little from ACDM model. How- 
ever, as mentioned earlier, the class of EU solution consid- 
ered here has provision for different composition of matter- 
energy in universe depending on the values of the parametes 
A and B. It can also accommodate a cosmological constant 
in a special case. The effective equation of state (We//) for 
EU remains negative always which is evident from fig. 6. 
The solid line in fig. 6 corresponds to the curve drawn using 
best fitted values. Dash and dotted curves are drawn with 
typical model parameters values within 95% and 99% confi- 
dence leves respectively. The transition of the universe from 
a deceleration phase to an accelerating phase in recent past 
is depicted from the curve of deceleration parameter against 
redshift plotted in fig. 7. The solid curve describes the one 
drawn with best fitted values and dotted and dash curves 
represent curves drawn with values within 99% and 95% 
confidence level respectively. We conclude that the class of 
EU solutions considered here is not ruled out by the obser- 
vations. However, this class of EU solutions admits different 
composition of matter-energies in the universe and the na- 
ture of composition depends on the value of parameter A 
in particular. The observations do in fact severely constrain 
the nature of allowed composition of matter-energy by con- 
straining the range of the values of the parameters for a 
physically viable model. 
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Figure 7. (Colour Online) Deceleration parameter (q) vs. redshift 
(z). Plot with the best fitted values of model parameters(Solid), 
with values within 95% confidence (Dashing) and within 99% 
(Dotted). 
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